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1. Introduction 

We prove |a;| _ ( d_2 ) decay for the critical two-point function of the continuous-time 
weakly self-avoiding walk in dimensions d > 4. This is a summary of the ideas 
and the steps in the proof. The details are provided in [12]. The proof is based 
on a rigorous renormalisation group argument. For the case d > 4, this provides 
an approach completely different from the lace expansion methods of |18l ITS] . But 
our main contribution is that our method applies also in the case of the critical 
dimension d = 4, where lace expansion methods do not apply. 

Renormalisation group methods have been applied previously to study weakly 
self-avoiding walk on a 4-dimensional hierarchical lattice. The continuous-time 
version of the model has been studied in the series of papers [H [T|)J [SJ H] ; see [5] for 
a review. More recently, a completely different renormalisation group approach to 
the discrete-time weakly self-avoiding walk on a 4-dimensional hierarchical lattice 
has been developed in [20] . 

The Ixl - ^ -2 - 1 decay for the two-point function for a continuum 4-dimensional 
Edwards model, with a smoothed delta function, has been proved in |24j ; unlike 
our model, this is not a model of walks taking nearest neighbour steps in the lattice, 
but it is expected to be in the same universality class as our model. The relation 
between our model and the Edwards model is discussed in [26| . A big step towards 
an understanding of the behaviour in dimension d = 4 — e is taken in |27| (their 
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model is formulated on a lattice in dimension 3 but it mimics the behaviour of the 
nearest-neighbour model in dimension 4 — e). 

Our renormalisation group method is a greatly extended and generalised form 
of work in [H|8l|9] for the hierarchical lattice and [T3l HH El [TT] for continuum quan- 
tum field theory. Details will appear in |12j . Our method is based on an exact 
functional integral representation of the two-point function of the continuous-time 
self-avoiding walk as the two-point function of a quantum field theory contain- 
ing both bosonic and fcrmionic fields. Such representations have been recently 
summarised in 1101. 



1.1. Background. A self-avoiding walk on the simple cubic lattice Z d is an 
injective map 

uj : {0,1,..., n} -> Z d (1) 

such that for all i, ui(i) and u>{i + 1) are nearest neighbours in Z d . We call n the 
number of steps. The main result of this article will actually be a statement about 
about random maps X : [0, T] — > Z d , but to explain the background we start with 
self-avoiding walk. 



Figure 1. An 8 step self-avoiding walk on Z d , d = 2. 



Let S n be the set of all self-avoiding walks with n steps and with ui(0) = 0. 
Let c n be the number of elements in S n . By declaring that all oj in S n have equal 
probability l/c„ we make S n into a probability space with expectation E„. The 
subscript n reminds us that the probability space depends on n. In the sequel 
"model" means a choice of probability space and law. 

This model arose in statistical mechanics. It is, for example, a natural model 
when one is interested in the conformation of linear polymer molecules. There is 
another natural model called the true or myopic self-avoiding walk. Unlike our 
model, true self-avoiding walk is a stochastic process which at each step looks at 
its neighbours and chooses uniformly from those visited least often in the past. 
Recent progress on this model is reported in [2"3"] . 

The key problem is to determine the growth in n of the mean-square displace- 
ment, 

E n \u,(n)\ 2 = c- 1 M™)! 2 ' ( 2 ) 



Self-avoiding walk in four dimensions 



3 



where |w(n)| is the Euclidean norm of u>(n) as an element of Z d . More precisely, 
we want to prove the existence of v such that 

lim n-' 2u E n \uj(n)\ 2 G (0,oo), (3) 

71— >OG 

and we want to calculate v. We will call this the v problem. 

As explained in [261 page 16], there is an easier version of this problem that we 
will call the Abelian v problem, because proving the existence of v after solving 
the Abelian problem is a Tauberian problem. Let S = lj n S n and let n(u>) = n for 
lo S n . For z > we define the two-point function 



uies 

Let 



G x (x) = 22z n Ml u(Mu))=m . (4) 



x (p) = J- z n M\ U (n{u))\* = J2 G ^)M P - ( 5 ) 

The Abelian version of the v problem is to determine the growth of y^( 2 ) /x^ as 
z j" z c , where z c is the common radius of convergence of the power series in this 
ratio. If v exists then it equals the Abelian v. In dimensions d > 5, according to 
the following theorem, v = 1/2. 

Theorem 1.1. [21, 22] For d > 5, there are positive constants A, D, c, e such that 

c n = A f i n [l + 0(n-% (6) 
E„| W (n)| 2 = £n[l + 0(n- e )], (7) 

and the reseated self-avoiding walk converges weakly to Brownian motion: 

w(LntJ) - B t . (8) 



VDn 

Also as \x\ — > oo, 

G Zc (x)=c\x\- (d - 2) [l + 0(\x\-% (9) 

The limit in © is called a scaling limit. The identification of scaling limits for 
dimensions d = 2,3,4 is the grand goal, but the v problem is a key intermediate 
objective because n~ u to(\nt\) is the candidate sequence for the scaling limit. 

If we set up the probability space without imposing the injective condition in 
the definition of uj, then the mean-square displacement is exactly n, because then 
the law for u is that of simple random walk. According to Donsker's Theorem, 
the scaling limit of simple random walk, with D = 1, is also Brownian motion. 
Thus, in dimensions d > 5 self-avoiding walk and simple random walk have the 
same scaling limit. When different models have the same scaling limit, we say 
the models are in the same universality class. One of the goals of mathematical 
statistical mechanics is to classify universality classes. 
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Theorcm ll.il will not hold with v =1/2 for dimensions four and less. There is 
a conjecture going back to [2] that, for d = 4, 



This and the next paragraph motivates our interest in four dimensions. 

In dimension d = 3, nothing is known rigorously about the v problem. The 
existence of v is not proved. It is not known that self-avoiding walk moves away 
from the origin faster than simple random walk, E„|u;(n)| 2 > n, nor is it known 
that self-avoiding walk is slower than ballistic, E n |w(n)/n| 2 — > 0. In dimension 
d = 2, there is the same basic lack of control as in d = 3, but the good news is that 
there is a candidate for the scaling limit, which tells us that if v exists it should 
be equal to 3/4. In [25], the process known as SLE 8 / 3 is identified as the scaling 
limit of self-avoiding walk subject to the unproven hypothesis that the scaling limit 
exists and is conformally invariant. 

SLE is a breakthrough discovery because it is provides a comprehensive list 
of possible scaling limits in d = 2. It has separated off the issues of existence of 
limits and universality and made it possible to study candidate limits without first 
having proved they are limits. On the other hand, theoretical physicists have a 
profound calculus called the Renormalisation Group (RG) that naturally explains 
when different models are in the same universality class and that can also prove 
the existence of limits. We will follow this path. RG, in the form that we will 
develop, was largely invented by Ken Wilson [28j [3Q1 ES]. RG as a rigorous tool 
originated with p] [15] . Later developments are reviewed in [6] . The hierarchical 
lattices mentioned earlier have special properties that greatly simplify RG. The 
n(logn) 1 / 4 growth of (fTUl) has been shown to hold for continuous-time weakly self- 
avoiding walk on a four dimensional hierarchical lattice in [4J |8j [9] . Very recently, 
the corresponding Abclian v problem has been solved in (20) for a discrete-time 
model on the hierarchical lattice. 

1.2. Continuous-time weakly self-avoiding walk and the main 
result. We now describe a probability law on a space of maps X : [0, T] — > Z d . 
We use the word "time" for the parameter t £ [0,T], but as for the discrete-time 
case there is a different space and law for each T. It is not a stochastic process 
which reveals more about itself as "time" advances, so it is better to think of the 
interval [0,T] as a continuous paramctrisation of a path in 1 d . 

Fix a dimension d > 4. Let X be the continuous-time simple random walk on 
Z rf with Exp(l) holding times and right-continuous sample paths. In other words, 
the walk takes its nearest neighbour steps at the events of a rate-1 Poisson process. 
Let P a and E a be the probability law and the expectation for this process started 
in X(0) = a. The local time at x up to time T is given by 



Afj, n {logn) 1/4 , E n \uj(n)\ 2 ~ £>n(log n) 1 / 4 . 



(10) 




(11) 



and we can measure the amount of self-intersection experienced by X up to time 
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1(0, T)= f dsi / ds 2 l x(si 
Jo Jo 



)=X(8 2 ) 



f f 

Jo Jo _^~ A 



2 

x,T- 



(12) 



Then, for g > 0, e 9 7 (°> T ) is our substitute for the indicator function supported on 
self-avoiding X. For g > 0, we define a new probability law 

P„, a (A) = E a (e-^ T H A )/E a (e-^ ^) (13) 

on measurable subsets A of the set of all maps X : [0,T] — > Z d with X(0) = a. 
For this model there is a ^ problerrp, but only the Abelian v problem for Z d is 
currently within the reach of the methods of this paper. 

The continuous-time weakly self-avoiding walk two-point function is defined by 

/>oo 

G g , v (a,b)= E a (e-^°< T h x{T)=b )e-» T dT, (14) 
Jo 

where v is a parameter (possibly negative) which is chosen in such a way that the 
integral converges. For p > define 

4 P V)= J2G g , u (a,b)\b-a\P. (15) 

bez d 

By subadditivity, c.f. |26) . there exists v c — v c {g) such that Xg\ v ) < oo if and 
only if v > v c . We call this v c the critical value of v. Our main result is the 
following theorem. 

Theorem 1.2. Let d > 4. There exists g m ax > such that for each g £ [0, c/ ma x] 
there exists c g > such that as \a — b\ — > oo, 

G 9 ,M 3 )M)= \ a \-2 + (16) 

This is the analogue of © in Theorem 1 1.1[ but now including dimension d = 4. 
There are no log corrections. Log corrections are only expected in the singular 
behaviour of Xg iy) as v \. v c for p > 0. The case g = is a standard fact about 
simple random walk; our proof is given for case g > 0. 



2. Finite volume approximation 

In this section we describe the first step in our proof, which is to approximate the 
infinite volume 1 d by finite volume, namely a discrete torus. 



1 solved on the hierarchical lattice for g small in [4j [8] [9] 
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We do not make explicit the dependence on g, which is fixed and positive. Let 
R > 3 be an integer, and let A = Z rf /i?Z d denote the discrete torus of side R. For 
a, b e A, let 

G A ,,(a, b) = J E a.A (e- 9J(0 ' T) lx(T)= b ) e~ vT dT, (17) 

where E a< \ denotes the continuous-time simple random walk on A, started from 
a. The following theorem shows that it is possible to study the critical two-point 
function in the double limit in which first A j" 7L d and then v \. v c . We will 
follow this route, focusing our analysis on the subcritical finite volume model with 
sufficient uniformity to take the limits. 

Theorem 2.1. Let d > 1 and v > v c . Then 

G v (a,b) = lim lim G Au ,(a,b). (18) 

v'iv A^Z d 



3. Integral representation 

The next step in the proof is to represent the two-point function in finite volume 
by an integral that we will approximate by a Gaussian integral. 

Recall that A denotes a discrete torus in 1 d . Given ip G C A and writing 
<p = i f x), i£A, we write dip x and d(p x for the differentials, we fix a choice of the 
square root \Z2ni, and we set 

Ipx = f— dkpx, i>x = —J=&Px- ( 19 ) 

Define the differential forms 

r x = tp x tPx + 4>x A V» x (x G A), (20) 

and 

ta,x = - ((Pit(- A0i + (-A^j;^ + ^ A (-A^ + (-A^) x A ^) , (21) 

where A is the lattice Laplacian on A defined by Ay^ = ^2 y .\ y - x \—i(<Pv ~ fx), 
and A is the standard wedge product. From now on, for differential forms u, v, we 
will abbreviate by writing uv = u A v. In particular ipxipy = — ''Py^Px and likewise 
tp x anticommutes with ip y and with ip y . The proof of the following proposition is 
given in [4j |9j ; see also [10| for a self-contained proof. 

Proposition 3.1. Given g > 0, let v be such that G\ tl/ (a,b) is finite. Then 

G\, v (a, b) = [ e -E. 6 A(TA 1 .+ S r. a +^r-)^ oy , 6 . ( 22 ) 
Jc A 
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The definition of an integral such as the right-hand side of (|2"2"|) is as follows: 

1. Expand the entire integrand in a power series about its degree-zero part 
(this is a finite sum due to the anti-commutativity of the wedge product, 
and the order of factors in the resulting products is immaterial due to the 
even degree), e.g., 

e - VT " = e -w.9.-£t*p.<*p. = e -"!P.e. ( i _ — dMfcdSx ) . (23) 

\ 2iri J 

In general, any function of the differentials is defined by its formal power 
series about its degree-zero part. 

2. Keep only terms with one factor ip x and one ip x for each i£A, write tp x = 
u x + iv x , fix = u x — iv x and similarly for the differentials. 

3. Rearrange the differentials to II^eA du x dv x , using the anti-commutativity of 
the wedge product. 

4. Finally, perform the Lebesgue integral over R 2 I A L 

This is explained in more detail in [10]. These integrals have the remarkable self- 
normalisation property that 

/ e - E .„ ( ..™,+^W., = 1 , ,24) 

Self-contained proofs of this, and of generalisations, can be found in [10] . The 
variables tp x and the forms ip x arc called fields. 



4. Quadratic or Gaussian approximation 

The integral representation of Proposition 13.11 opens a natural route for approxi- 
mation by non-interacting walk with different parameters. To do this we split the 
exponent ta.x + 9T X + vt x in (|22[) into a part which is quadratic in the variables ip 
and a remainder. When the remainder is ignored the rest of the integral becomes 
Gaussian and the Gaussian integral represents a non-interacting walk. It is impor- 
tant not to assume that the best approximation is the quadratic terms ta.x + vt x . 
We even want to allow t\ to be divided up. To see what a different coefficient in 
front of ta means we make the change of variable tp x i— > s/1 + ZQtp x , with zq > — 1. 
This gives 

Ga,„M) = (1 + *o) / e-^^( (1+ ^ )rA ^^ 9(1+2o)2T " +iy(1+Zo) ^)^ a ^, (25) 
Jc A 

where the Jacobian is contained in the transformation of ip^ip. Then, for any 
m 2 > 0, simple algebra allows us to rewrite this as 



G A ,„M) = (l + Zo) / e- S(A) ^ o(A W b , 



(26) 
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where 

S(A) = ^(r A ,* + mV), (27) 
xeA 

(A) = (do^ + "0T X + Z T A , X ) , (28) 



go = (1 + zq) g, i/ = (1 + z Q )v c , m = (1 + zq)(v - v c ), (29) 

and v c was denned below (|15j). The two-point function G^ v (a, b) in (|26|) does not 
depend on (zo,«i 2 ) so, in the next theorem, these are free parameters that do not 
get fixed until Section [TJ1 In view of Theorem 12.11 and Proposition 13.11 to prove 
Theorem 11.21 it suffices to prove the following theorem. 



Theorem 4.1. Let d > 4. There exists c/ m ax > such that for each g £ [0, <? m ax] 
there exist c{g) > such that as \a — b\ — > 00, 

lim lim (1 + zo) [ e- 5 ^-^ A W& = i ^1 2 (1 + o(l)) . (30) 
i4"cAfZ d J c a \a-b\ d - ly 



To prove Theorem 14. 1[ we study the integral on the left-hand side via a renor- 
malisation group analysis, without making further direct reference to its connection 
with self-avoiding walks. In order to calculate this integral we define, for a G C, 

Vb(A) = V (A) + atp a + atp b (31) 

and use 

^sw-vow _ dd I e _ s(A) _ MA) ^ (32) 

00 da 



Jc A 



We will call a an external field. 



5. Forms and test functions 

In this section we introduce notation for handling the differential forms that appear 
in Theorem 14. II We will write form in place of "differential forms" from now on. 
We focus on dimension d = 4, but leave d in various formulas since 4 can also 
appear for other reasons. 

5.1. The space A/". A form is a polynomial in ip,^ with coefficients that are 
functions of (if, a) e C A x C. 

Given a £ C we define = a and ct 2 = ^ so that a can be identified with a 
function er : {1,2} — > C. Similarly, let A 2 = A x {1, 2} so that given if £ C A we 
have the function on x = (s, i) £ A 2 defined by 



(33) 
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Since </> and ip are in one to one correspondence and since we are only interested 
in functions on A2 that arise from some <p we write (j) G C A . 

Forms are elements of the algebra Af whose generators are the degree one forms 
(ipxi'tpx, x € A) subject to the relations that all generators mutually anticommute. 
For x = (s, i) G A2, we write 

*■-{* !:i 

Then we introduce the space A* = U^ A| of all sequences in A 2 with finitely 
many terms so that every monomial in ip can be written in the form, for some 

y g A*, 

{Vyi ■ ■ ■ Vy, if?>l- 

The q = term in A* is a set consisting of a single element called the "empty 
sequence", which by definition has length zero. Given a sequence y G A*, q = q(y) 
is the length of the sequence and yl = q{y)\- Every element of Af has the form 

F = F(<f>,a) = ^FyfaW. (36) 
Given x — (xi, . . . , x p ) G A 2 and z — (z\,... ,z r ) G {l,2} r , we write 

For X C A, we define Af(X), which is a subspace of Af, by 

Af(X) = {F G Af : F x>y = if any component of x, ?/ is not in X}. (38) 

For example t x G 7V({x}) and ta, x G Af(X) where X = {y : \y — x\ < 1}. 
By introducing 

we write the formal Taylor expansion of F ((f) + £) in powers of £ and a as 

K,yeA*,ze{l,2}* y 

Functions / : A* x A* x {1, 2}* — > C are called iesi functions. We define a pairing 
between elements of Af and the set of test functions as follows: for a test function 
/, for $ G C A , let 

(fju= E tAtt/-.-"-^./:...-- (41) 

a:,?/gA*,ze{l,2}* 



(39) 
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5.2. Local polynomials and localisation. For a function / : A -> C 
and e a unit vector in Z d we define the finite difference derivative (V e /) K = 
/(a; + e) — f(x). Repeated differences such as (V e V e '/) 2 . are called derivatives. 

A local monomial is a product of finitely many fields and derivatives of fields 
such as M = ipipV e f. Using this example to introduce a general notation, given 
x £ A let M x = ipx'fpxC^e^x, and given X C A let M(X) = J2xex^x- Local 
polynomials are finite sums of local monomials with constant coefficients. 

An important example of a local polynomial is 

V = gr 2 + vt + zt^ x + Xl a af + Xl^a(p + qaa, (42) 

which extends the local polynomial of (f3"Tj) by the addition of the term qaa. The 
indicator function 1 Q : A — > {0, 1} equals 1 when evaluated on a and is zero 
otherwise. The parameters (g, v, z, A, q) are called coupling constants. 

Euclidean symmetry: The lattice Z d has automorphisms E : 7h d — > Z d . An 
example for d = 1 is Ex = 1 — x. By letting an automorphism E act on the 
spatial labels on fields, ip x i— > ifEx, E induces an action, E : Af — >• A/". A local 
polynomial P is Euclidean invariant if automorphisms of Z d that fix x also fix . 
For example, ipipVe'P is not Euclidean invariant because there is a reflection that 
changes (p x + e into tp x -e so that (V e (p)x >-> (V-e^)^. On the other hand, the term 
ta in (|42[) is a Euclidean invariant local monomial. 

Gauge invariance: A local polynomial is gauge invariant if it is invariant under 
the gauge flow: (cr,ip) — > (e ie a,e l6 tp). Thus V of (|4*2"|) is gauge invariant. 

Super symmetry: There is an antiderivation Q : J\f — > N characterised by 

Q^x = ipx, Qi>x = - fx, Qfx = ipx, Q^Px = fx- (43) 

An element of F £ Af is said to be super symmetric if QF = 0. The terms r, ta, t 2 
in V are supcrsymmetric local monomials. The forms af, af, aa are gauge invari- 
ant, but not supcrsymmetric. It is straightforward to check that Q 2 generates the 
gauge flow. Therefore supcrsymmctry implies gauge invariance. Further details 
can be found in [ID] . 

The pairing (|4T|) defines F £ Af as a linear function, / H> (F,f) , on test 
functions. The subscript means that we set = 0. Let II be a set of test functions. 
Two elements F\ and F2 of Af are equivalent when they define the same linear 
function on II. We say they are separated if they are not equivalent. 

Example 1. Let II be the set of test functions that are linear in their A arguments. 
Fix a point k e Z d . Let F = fk, and let F' = fo + (k ■ Vf)o- Then F and F' are 
equivalent because a linear test function f(x) = a + b ■ x cannot separate them: 

(FJ)=a + b-k, (F'J) =a + k-Vf = a + b-k. (44) 

To avoid confusion let us emphasise that two different contexts for "polynomial" 
are in use: a test functions can be a polynomial in x £ A, while local polynomials 
are polynomial in fields. 
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The choice for II in this example is not the one we want. The details in the 
definition given below are less important than the objective of the definition, which 
is that II should be a minimal space of test functions that separates the terms in 
621). 

We define II to be the set of test functions f(x,y,z) that are polynomial in 
the A arguments of (x, y) £ A* x A* with restrictions on degree listed below. For 
/ £ II, as a polynomial in the x, y components in A: 

1. The restriction of / to (x,y,z) with r(z) = has total degree at most d — 
p(x)[4>] - q{y)[(f>}; f(x,y,z) = when d - p{x) [<f>] - q{y)[<j>\ < 0. Here 

[0] = (d-2)/2. (45) 

For dimension d = 4, [0] = 1. 

2. The restriction of / to (x,y,z) with r(z) = r G {1,2} has total degree at 
most r — p{x) — q(y); f(x, y, z) — if r — p{x) — q(y) < or r > 2. 

Let V be the vector space of gauge invariant local polynomials that are separated 
by II and, for X C A, let V{X) = {P{X) : P G V}. The following proposition 
associates to any form F € N an equivalent local polynomial in V(X) [12] . 

Proposition 5.1. For X C 7L d there exists a linear map Locx '■ M — > V(X) such 
that 

(a) (L^ x F,f)o = (FJ) for f e II, F e Af, (46) 
(6) F(Loc x F) = Imc E x(EF) for automorphisms E : Z d -4 if, F G TV, (47) 
(c) Locx- o Locx = Locx' /or X, X' C A. (48) 

Let Vh C V be the subspace generated by monomials that are not divisible by 
a or <7, and let Vo C V be the subspace generated by monomials that are divisible 
by a or a. Then V = Vh © Vo, and on this direct sum we define 

Loc x = Loc x © Locxn{o,6}- (49) 

Symmetry considerations for the integral representation restrict the domain of Loc 
in our applications so that its range reduces to polynomials of the form V as in 
62). 



6. Gaussian integration 

6.1. The super-expectation. For a A x A matrix A, we define 

Sa(A) = ^ (vxAcy&x + *l>xAxv$v) - ( 50 ) 
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When A = m 2 — A this is the same as S(A) which was defined in (p?T|) . Let C be a 
positive-definite A x A matrix. Then A = C^ 1 exists. We introduce the notation 

E C F = { e- SA(A) F, (51) 

for F a form in Af. The integral is defined as described under Proposition 13. II We 
call C the covariance because Ec<^ a 06 = C a b- More generally, if F is a form of 
degree zero, i.e., a function of <f>, then EpF is a standard Gaussian expectation for 
a complex valued random variable <p with covariance C |10j . 

We define a space Af x in the same way as Af is defined, but with <f> doubled 
to (4>,£) so that (4>,?P) doubles to the pair (0, VOi (^i 7 ?) with rj = {2-Ki)~ 1 / 2 d^. The 
external field a is not doubled. We define : N — >• A/" x by 

(0F)(<M) = E 3W + 0(</> + *?)"■ (52) 

yeA« V ' 

We write Ec#-F for the element of Af obtained when the integral over C A in Ec 
applies only to (£,,r/). In the general case where F is a form this is not standard 
probability theory, because McOF takes values in Af. To keep this in mind we call 
this a super- expectation. The variables and forms (£, 77) that are integrated out are 
called fluctuation fields. 

6.2. Finite-range decomposition of covariance. Suppose C and 
Cj , j = 1 , . . . , N, are positive-definite A x A matrices such that 

C = Y,Cj. (53) 
3=1 

Let C = £fe=2 C fc . Then, as in the standard theory of Gaussian random variables, 
the Ec expectation can be performed progressively: 

E C F = E c >+ Cl F = E C ' (E Cl OF) . (54) 

For further details, see [12]. 

From now on we work with C = (m 2 — A) -1 , where A is the finite difference 
Laplacian on the periodic lattice A. Given any sufficiently large dyadic integer L, 
there exists a decomposition C = Ylj=i Cj such that Cj is positive-definite and 

Cj(x,y) = if \x-y\>V. (55) 

This is called the finite range property. The existence of such a decomposition is 
established in [7] for the case where A is replaced by 7L d . In [SJ Lecture 2] it is 
briefly explained how the decomposition for the periodic A case is obtained from 
the Z d case, for A a torus of side L N . To accommodate this restriction on the side 
of A the infinite volume limit in Theorem 14. II is taken with a sequence of tori with 
sides L N , N e N. 
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We conclude this section with an informal discussion of scaling estimates that 
guide the proof. Equation (|54|) says that F, which depends on a field with covari- 
ance C, can be replaced by E^OF, which depends on a field characterised by the 
covariancc C . Repeating this operation j times will replace F by a new F that de- 
pends on a field at scale j characterised by the covariance X)fcLj+i According 
to estimates in [7] , this sum is dominated by the first term which satisfies 

\^ c ^yG ]+ i{x,y)\ < const i-VM-Mii-I^W, (56) 

where the symbol [</>], which is called the dimension of the field, was defined in 
(|45p. The typical field at scale j behaves like "half a covariance," and in particular 
the standard deviation of <p x is w L~^'. Furthermore, the estimate on derivatives 
in (|56p says that typical fields at scale j are roughly constant over distances of 
order LP . 

We can now explain why the terms in V as defined by (|42[) play a pre-eminent 
role. For a cube B of side U , which contains L dj points, 

J2 <$, x « L^ d -P^. (57) 

In the case of d = 4, for which [</>] = 1, this scales down when p > 4 and <p p is said 
to be irrelevant. The power p = 4 neither decays nor grows, and is called marginal. 
Powers p < 4 grow with the scale, and are called relevant. Since the derivatives in 
(|56|) provide powers of L, the monomial ip{— A)<p is marginal. Thus r, ta,t 2 are 
the supcrsymmetric marginal and relevant monomials. 

6.3. Progressive integration. To prove Theorem |4~T1 using ([32]) we have 
to calculate 

f e -S(A)-V (A) =Ece -V (A) j (58) 

where Vq is given by (|3T|) . This Vq equals V as defined in (|4"2|) . with (g, v, z, A, q) 
replaced by (g , i/ , z , A , go) with 

go = 0, A = 1. (59) 

Sections 16.11 and 16.21 have taught us that we can evaluate Ece" Vo ' A ' by the 
following iteration: let 

Z a = e- Vo(A) . (60) 
Inductively define Z } ;, j = 0, . . . , N, by 

Z J+1 =E Cj+1 ^, (61) 

Then 

E c e- y ° (A) = Z N . (62) 

Therefore the proof of Theorem 14.11 now depends on the analysis of the sequence 
Zj. Our proof will depend on showing that the Zj simplify as j increases. In fact, 
in the next section we will see that they become more Gaussian, in the sense that 
the gr 2 term becomes smaller. The index j will be called a scale. 
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7. Perturbation theory and flow equations 

In this section we start to prove that Zj becomes more Gaussian as j increases. To 
do this we adapt to our particular setting a pcrturbative calculation of the kind 
that appears in [30] . 

For X C A and V as defined in (|32"]l. define 



where 



with 



I j MV)=e- v ™{l + ±W j (V,X)), (63) 
W 3 (V,X) = (1 - Loc x )F Wj (V(X),V(A)) (64) 

3 

Wj =Y,Ci, (65) 
F W] (V(X),V(A)) = J2 ^(D n R V(X))w?(DlV(A)); (66) 

n>l 

the latter sum truncates at n — 4 due to our quartic interaction. The symbols Dr 
and Dl denotes right and left differentiation with respect to fields. The "left/right" 
is to specify signs, but this and the precise definition arc not of immediate impor- 
tance, so we just give an example. 

Example 2. For V = 44 and X = {x}, (D%V(X))w] (Z>£^(A)) equals 

\-T, y eA w jfav) n = 2 - 

When j = 0, Ij,x(V) — e~ v ^ because wq — 0. Therefore we can choose 
the coupling constants to make it equal to Zq. Furthermore, Ij,x{Y) has the 
martingalc-likc property exhibited in Proposition [7T] which says that integrating 
out the fluctuation field is approximately the same as changing the coupling 
constants in V to new coupling constants called (<? pt , v pt , z pt , A pt , q P t)- The formu- 
las for the new coupling constants are called perturbative flow equations. 

Proposition 7.1. As a formal power series in (g, v, z, A, q), 

Ec j+ J j>A (V) = I j+ i,A(Vpt) mod (g,u,z,X,q) 3 , (68) 



where 



V pt = V pt (V) 



(69) 
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has the same form (|42p as V, with (g,i>,z,\,q) replaced by 

9 P t = 9 - c g g 2 + r pt v (70) 



i/^ = i/ + 2^ i+1 (0,0)+r£., (71) 

I pt 



z + rfj, (72) 



l+gS j + -(ML 

q pt = q + \ 2 C J + 1 (a,b), (74) 

where c g > and 

AT 

M = L 2 ^ + 2 9 ^ C fe (0,0)), (75) 
fc=i+i 

and where r^* • , • and j are computable uniformly bounded homogeneous poly- 
nomials of degree 2 in (g, fi, z). There are g 2 terms in j, but they are summable 
in j and therefore do not overpower c g g 2 ; Sj is a summable sequence of positive 
numbers determined by Cj . 



The (3 function. The right hand side of ((70)) is known as the f3 function. The 
simpler recursion obtained by setting = 0, namely 

ffi+i = 9j ~ Cgg 2 , g = go, (76) 

creates a sequence g~j that tends to zero like j _1 as j — } oo. The sequence Zj 
becomes more Gaussian due to the famous observation, known as infra-red asymp- 
totic freedom, that (|76p controls the behaviour of the more complex recursion of 
Proposition 17. II and drives the r 2 term to zero. 



8. The renormalisation group map 

The problem with the second order perturbative calculation in Section [7] is that 
the error is not only of order 3 in the coupling constants. It is also at least of 
order 3 in L~ d: >\A\. In fact there will also be an exp(0(L~ dj \ A|)) in this error! 
We will call this the exp(0(|AQ) problem. The remedy is not to work with Ij(A), 
but with IIbca where B is a cube and the allowed cubes pave A. The idea 
is that by choosing the side of B to be bigger than the range of Cj+i, we can 
take advantage of independence of cubes that do not touch to more or less use 
our perturbation theory with A replaced by individual cubes. This idea requires a 
systematic organisation which we describe in this section. 
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8.1. Scales and the circle product. Let L > 3 be an integer. Let 
R = L N , and let A = 1 d /{R1 d ). 

Definition 1. (a) Blocks. For each j = 0,1, . . . , N , the torus A is paved in a 
natural way by L N ~i disjoint d-dimensional cubes of side LP . The cube that 
contains the origin has the form (for L odd) 

{x G A : \x\ < hp 1)}, (77) 

and all the other cubes are translates of this one by vectors in ]Jl, d . We call these 
cubes j-blocks, or blocks for short, and denote the set of j-blocks by Bj = Bj(A). 

(b) Polymers. A union of j-blocks is called a polymer or j-polymer, and the set 
of j-polymers is denoted Vj = Vj(A). The size \X\j of X G Vj is the number of 
j-blocks in X. 

(c) Connectivity. A subset X C A is said to be connected if for any two points 
x a ,xi, G X there exists a path (xi,i = 0, 1, . . . n) G X with — ajj||oc = li = 
x a and .t„ = Xb- According to this definition, a polymer can be decomposed into 
connected components; we write C(X) for the set of connected components of X. 
We say that two polymers X, Y do not touch if min{||a; — j/||oo : x € X, y e Y} > 1. 

(d) Small sets. A polymer X G 'Pj is said to be a smaH sei if < 2 d and X is 
connected. Let Sj be the set of all small sets in Vj. 

(e) Small set neighbourhood. For X C A let 

A*= |J K (78) 

The polymers of Definition [T] have nothing to do with long chain molecules. 
This concept has a long history in statistical mechanics going back to the important 
paper [17] . 

Proposition 8.1. Suppose that X-y, . . . ,X n G Vj+i do not touch each other and 
let Fi(Xi) G Af(Xi). The expectation Ecj+i has the factorisation property: 

n n 

Ec 3+1 HFi(Xi) = [] Vc i+1 Fi(Xi). (79) 

i— 1 m— 1 

Proof. Gaussian random variables are independent if and only if the off-diagonal 
part of their covariance matrix vanishes. This generalises to our forms setting, and 
so the proposition follows from the finite range property of Cj+i. □ 

Given forms F, G defined on Vj , let 

(FoG)(A)= F(X)G(A\X). (80) 

xevj 

This defines an associative product, which is also commutative provided F and G 
both have even degree. 
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8.2. The renormalisation group map. Recall that we have defined 
Ij,x(V) in (|63|) . Given a yet-to-be-constructed sequence Vj, for X £ Vj, let 

i 3 {x)= n ijAVj)- (si) 

BeBj 

We have defined V in ((3TJ). Let K (X) = l x =z- Then the Z defined in (JBDJ) is 
also given by 

Z = I (A) = (I o aK )(A), (82) 
because Jo, A (A) = e~ y °' A ) since wo = 0. 

Definition 2. We say that A' : Vj AT has the component factorisation property 
if 

A-(X)= J] A(F). (83) 
Yec(x) 

Suppose, inductively, that we have constructed (Vj, Kj) where Kj :Vj^N is 
such that 

(<) Zj = (Ij oKj)(A), 

(ii) Kj has the component factorisation property, (84) 
{Hi) For X £ Vj, Kj(X) £ Af(X*). 

Our objective is to define (Vj+i, Kj + i), where Kj + i : Vj+i — > Af has the same 
properties at scale j + 1. Then the action of l&c j+1 on Zj has been expressed as 
the map: 

(Vj,Kj)^(V j+1 ,K j+1 ). (85) 

This map will be constructed next. We call it the renormalisation group map. 
Unlike Zj <— > ~E8Zj it is not linear, so this looks like a poor trade, but in fact it is a 
good trade because the data (Vj, Kj) is local, unlike creatures such as exp(— Vj(A)) 
in Zj. The component factorisation property and Proposition l8 . 1 1 allows us to work 
with Kj on the domain of all connected sets in Vj. We can prove that Kj(X) is 
very small when the number of blocks in X is large; in fact, only the restriction of 
Kj to the small sets Sj plays an important role. 



9. The inductive step: construction of Vj+± 

In accordance with the program set out in Section 18.21 we describe how Vj+i is 
constructed, given (Vj,Kj). Our definition of Vj+i will be shown to have an 
additional property that there is an associated Kj+i, which, as a function of Kj, 
is contractive in norms described in Section [TUl 
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Recall that the set S of small sets was given in Definition [T] For B 6 Bj define 
Vj+i to be the local interaction determined by: 

V 3 {B) = Vj{B) + Loc B ]T ]V-M Y r lK AY)> 

YCS^YDB 1 * 1 (86) 

where V pt = V pt {V) with generic argument V is defined in (f69|). Recalling the 
discussion of "relevant terms" just after (j57|) . one sees in (|86|) that Vj+i has been 
defined so that relevant terms that would expand if they remained inside Kj are 
being absorbed into Vj+i. 

We have completed the V part of the inductive construction of the sequence 
(Vj,Kj). Before discussing the K induction we have to define some norms so that 
we can state the contractive property. 



10. Norms for K 

Let f)j > and Sj > 0. For a test function / as defined in Section HT2"1 we introduce 
a norm 

11/11*,= sup sup f)J p - q sJ r L^\V a f x , y J. (87) 

x,j/eA*,2e{i.2}* HooO 

Multiple derivatives up to order 3 on each argument are specified by the multi-index 
a. The gradient V represents the finite-difference gradient, and the supremum 
is taken componentwise over both the forward and backward gradients. A test 
function / is required to have the property that f x ,y,z = whenever the sequence 
x has length p > 9 or the sequence z has length r > 2; there is no restriction on 
the length of y. By the definition of the norm, test functions satisfy 

\V a f x ,y,*\ < bf^-^i/ll*,- (88) 

We discuss the choice of Sj in Section [12] when it first plays a role, and here we 
focus on f)j. An important choice is 

^=l j= g a L-^\ (89) 

for a given £q- The is there because unit norm test functions of one variable 

should obey the same estimates as a typical field, and test functions of more than 
one variable should obey the estimates that a product of typical fields obeys. 
Recall the pairing defined in (jlTj) and, for F £ J\f and (f> 6 C A , let 

11^11^.= sup \(F,g)t\. (90) 

9 : llsll<fj <! 

The following proposition provides properties of this seminorm that are well adapted 
to the control of K. 
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Proposition 10.1. Let F,Fx,Fz G M. The norm obeys the product property 

ii^i^Ht^ < mw^miT^, (91) 

and, if £q is chosen large enough, the integration property 

mc i+ Mn AW <^c j+ A\F\\T, + , A 2 W - (92) 

For further details, see [12]. The second conclusion shows that the norm controls 
the forms when a fluctuation field is integrated out: on the right hand side the 
norm is a zero degree form, and hence the expectation is a standard Gaussian 
expectation. 

The most important case of the seminorm is the case <f> = 0, but knowing 
that (X)||t < oo cannot tell us whether K(X) is integrable. For this we must 
limit the growth of K(X) as <j) ~ > 00 > arL d the resolution of this issue will be 
obtained using Definition [3] below. 

Our intuitive picture of Kj(X), where X G Vj, is that it is dominated by a 
local version of the remainder (g, v, z, A, q) 3 in (|68|) . To estimate such remainders 
we must, in particular, estimate Ij : x which contains exp(— g 3 - ^2 x£X IVx| 4 )- By 
(|56|) the typical field ip at scale j is roughly constant on scale U , and X contains 
0(U d ) points. Therefore this factor looks like cxp(— gjL dj |<p| 4 ). This is a function 
of Lp/\)j with f)j w g. 1 / 4 i _ J d / 4 j which in four dimensions can be rewritten as 

g- 1 / 4 L _J ^] because [(/)] = 1. We want to prove that gj decays in the same way 
as does g~j in (|76|) . and with this in mind we replace gj by the known sequence gj. 
This leads us to our second choice 

t),=^=fc o .g7 1/4 ^ [01 , 

where the constant kg is determined so that exp(— Vj(B)) will, uniformly in j, have 
a T^(hj) norm close to one. 

In the previous discussion we made the assumption that the typical tp at scale 
j is roughly constant on scales U> . Our norm recognises this; it is a weighted 
norm, where the weight permits growth as fields become atypical. The weight is 
called a large field regulator and is defined next. 

Consider a test function / that is an ersatz field (p, namely a complex-valued 
function / = f x for x € A. For X C A, we write / € n(AT) if / restricted to X is 
a polynomial of degree three or less. We define a seminorm on <j) = (tp, (p) by 

UU ](X) = infill - fh jM : / G n(X)}; (93) 
note that we are setting f)j = £j in the above equation. 

Definition 3. Let j <G No, X <G Vj, and 4> G C A . The large-field regulator is given 
by 

Gj(X,cf>)= J] expM!. (B , ); (94) 
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where B* is the small set neighbourhood of B defined in ([75| . For each X G 
Vj, we define a seminorm on Af(X*) as follows. For K{X) G Af(X*), we define 
ll-K^X)!!^ h . to be the best constant C in 

\\K(X)\\ nAhj) < CGj(X,$), (95) 

where we have made explicit in the notation the fact that the norm on the left 
hand side is based on the choice f)j = hj. 



11. The inductive step completed: existence of Kj + i 

We have already specified Vj+i in Now we complete the inductive step by 

constructing Kj + i such that holds. The following theorem is at the heart 
of our method |12j . It provides Kj + i and says that we can continue to pro- 
long the sequence (Vj, Kj) for as long as the coupling constants (gj, i/j, zf) remain 
small. Moreover, in this prolongation, the To norm of -fiTy+i remains third or- 
der in the coupling constants and is therefore much smaller than the perturbativc 
(K- independent) part of Vj+i. 

For a > 0, set fj(a, 0) — 0, and define 

fj(a,X) = 3 + a(\X\j - 2 d ) + , X G Vj with X ^ 0. (96) 

Note that fj(a,X) = 3 when X G Sj, but that fj(a,X) is larger than 3 and 
increases with the size of \X\j if X^Sj. We fix a to have a sufficiently small 
positive value. 

The following theorem is proved for two different choices of the norm pairs 
|| • || j and || • || j+i, in (|97[) and (|98|) . and for two corresponding choices of the small 
parameter esi, as follows: 

• II • II, = II ' ll<5 Jlfc , with hj = k a g- 1/A L~^\ and || ■ = || ■ \\g j+1 ,h j+1 with 
hj+i = kog~j+l 4 . The small parameter esi is proportional to g^ 4 - 

• II ' \\j = II • IItoA- with *j = ZqL-1W, and || ■ || i+1 = || • ||t / j+1 . The small 
parameter esi is proportional to gj . 

Define a cone C = {(gj, Vj, Zj)\g > 0, |z^| V \z\ < bg, gj < c(b,L)}. The constant b 
is arbitrary, and c(b, L) is a function of b, L constructed in the proof of the next 
theorem. 

Theorem 11.1. Let (gj,Vj,Zj) G C. Let a be sufficiently small, and let M be 
any (large) positive constant that is independent of d,L. There is a constant c p t 
( depending on d,L) such that the following holds. Suppose that Kj : Vj — > Nj has 
properties (|84p and satisfies 

\\Kj(X)\\j < Mc pt e f s } {a ' X \ X G Vj connected, (97) 
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Then, if L is sufficiently large (depending on M), there exists Kj+i : Vj+i — > A/j+i 
with properties (|84|) at scale j + 1 and 

\\K j+1 (U)\\ j+1 < 2c pt e£ +l(a ' C/) , (7 G V j+1 connected. (98) 



12. Decay of the two-point function 

Finally, we combine the machinery we have developed, to outline the proof of The- 
orem [O] As we have already noted, Theorem ll.2l is a consequence of Theorem 14. II 

We must study the coupling constant flow. The linear map Loc b '■ A/" — > V is 
bounded in To norm |12j , so according to the inductive assumption (|97|) on the To 
norm of Kj, the coupling constants in Vj of ([86]) are small (third order) adjust- 
ments to the coupling constants in Vj. Theorem 1 1 1 . 1 1 ensures that this smallness 
is preserved as the scale advances. 

We first consider the case (Ao,<7o) = (0,0). In this case, (Xj,qj) = (0,0) for all 
j. The definition of Vj+i in (|86|) then gives rise to a non-perturbative version of 
the flow equations of Proposition 17. 1[ in which the effect of K is now taken into 
account. When Vj M> Vj+i is expressed as 

(gj,Uj,Zj) h-> (g j+l ,Vj +1 ,Zj +l ) (99) 

we find that 

9j+i = 9j - Cgffj + r g ,j , (100) 

v j+1 = Vj + 2gC j+x {Q, 0) + r vJ , (101) 

zj+i = Zj + r Ztj , (102) 

K J+ i = r K ,j(gj,Vj,Zj,Kj), (103) 

where the r's now depend also on Kj, and where we have added the map rx.j ■ 
(gj, Vj, Zj, K 3 ) i y Kj + i defined by Theorem lll.il Furthermore, we prove that the 
r's are Lipschitz functions of (gj, Vj, Zj, Kj), where K belongs to a Banach space 
normed by a combination of the norms in Section 111! These are the properties 
needed to prove that K only causes a small deformation of the perturbative flow 
V^V pt . 

The main theorem now reduces to an exercise in dynamical systems. We prove 
that there is a Lipschitz stable manifold of initial conditions (zq,vq) = h(m 2 ,go) 
for which the sequence (Vj, Kj), j — 0, . . . , N , has a limit as N — > oo and m 2 I 0. 
We call this the global trajectory. For m? = 0, the global trajectory tends to the 
fixed point (V,K) = (0,0). In particular, gj — > 0, which is infra-red asymptotic 
freedom. Referring to (|29p . we have four unknown parameters gQ,VQ, zq, m 2 related 
by three equations, and now there is a fourth equation (zq,vq) = h(m 2 ,go). By 
the implicit function theorem we solve for the unknowns as functions of (g, v). As 
v I v c (g), m 2 4- and vice- versa. 
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Now we consider the flow for (\j,qj). According to (|55|) . A = 1 and qo = 0. 
Using (|49|) . we prove that the terms r g j,r u j, r z j do not depend on \j,qj and thus 
the coupling constants g, v, z have no dependence on A, q. From we find 



X j+ i =(l + 9i Si + \^L-^C { ^ A, + r AJ , (104) 

qj+i = qj + A 2 Cj+i (a, b) + r g>j , (105) 

where r\ t j,r q j are corrections arising from Kj. 

Recall that Sj was defined in Definition [TJ Let s a ,b be the first scale j such 
that there exists a polymer in Sj that contains {a,b}. The correction r q j is zero 
for all scales j < s a j,: according to (l49l) and the definition of V in (|86[) there can 
be no tTtJ contribution from Kj until the first scale where there is a set X g 6>j 
that covers {a, 6}. Also, by the finite range property, Cj+\ (a, b) = for j < s a j,- 
Thus (fTU5|) gives 

N 

1N= J2 (^(': ■ >;,,)■ (106) 

J=S<J,i, 

At scale A, A is a single block in Bn, so by the definition of the circle product, 
Zjq is simply given by 

Z N = {I N o K N )(A) = I N (A) + K N (A). (107) 

The final renormalisation group map is the action of Ecni not Mc N 0. This means 
that the fields </>, ij) are to be set to zero in In, Km, and only dependence on a 
remains. By (|63j) we compute two o derivatives of In and find 



if- 



dodo 



Z N = q N - K^, where A ff(T = — 
oerocr 



(108) 



The oo derivative is a coefficient in the pairing (|4"Tj) , and the T norm bounds this 
pairing, so Theorem 111.11 gives 

\K 9a \ < \\K\\ To s N 2 < 0(g? )s N 2 . (109) 

We are able to prove Theorem 111.11 with 

*j=*^L*, h ™0{L*">*), (110) 

where So is a constant, so that, when N > s aj b, 

\K Sa \ < 0(g%)L- 2N ^ = 0(g N \a - b\~ 2 ). (Ill) 

This tends to zero as N — > oo. 

By a similar estimate we can control the rxj, r q j terms in f|104[) . ()106[) . These 
contain o derivatives of the Kj terms in (|86j). The conclusion is that Aoo = 
limAr_ ! . 00 Atv and qoa = liniAr_ J . 00 qjy exist and are bounded away from zero. 
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By ([32]), the left hand side of flU is given by 

lim(l + zq) lim / e- s{K) - Vo{K) ^ b = ( lim(l + z )) lim 9oo . (112) 
From (fT04]) and (fT06|) wc find that 

OO 

lim^-A^ V C j+ i(a,b), (113) 

where m 2 = in Cj+i, and ~ means that the ratio of the left hand side and the 
right hand side tends to one as a — b — > oo. Next, we use the finite range property 
to restore the scales j < Saj, to the sum, which then becomes the complete finite 
range decomposition for the infinite volume simple random walk two-point function 
(-A)~\a,b), 

lim qoo ~A^(- A)- 1 (a,b). (114) 

m 2 \.0 

The right hand side of Q114p . and hence of (|112[) , is thus asymptotic to a multiple of 
\a — b\~ 2 as \a — b\ — > oo, as desired, since the inverse Laplacian has this behaviour. 
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